https://filesoo.com/product/8816/geometric-numerical-integration-structure-preserving-algorithms-for-ordinary-differential-equations